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Abstract. Belgium is amongst few artificial countries, established on pur- 
pose, when Dutch and French speaking parts were joined in a single unit. This 
makes Belgium a particularly interesting testbed for studying bio-inspired tech- 
niques for simulation and analysis of vehicular transport networks. We imitate 
growth and formation of a transport network between major urban areas in 
Belgium using the acellular slime mould Physarum polycephalum. We repre- 
sent the urban areas with the sources of nutrients. The slime mould spans 
the sources of nutrients with a network of protoplasmic tubes. The protoplas- 
mic tubes represent the motorways. In an experimental laboratory analysis 
we compare the motorway network approximated by P. polycephalum and the 
man-made motorway network of Belgium. We evaluate the efficiency of the 
slime mould network and the motorway network using proximity graphs. 
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1. Introduction 

Plasmodium is a vegetative stage of the acellular slime mould Physarum poly- 
cephalum. This is a single cell with many nuclei. The plasmodium feeds on mi- 
croscopic particles |40j . During its foraging behaviour the plasmodium spans scat- 
tered sources of nutrients with a network of protoplasmic tubes. The protoplasmic 
network is optimised to cover all sources of food and to provide a robust and 
speedy transportation of nutrients and metabolites in the plasmodium body. The 
Plasmodium's foraging behaviour can be interpreted as computation. Data are 
represented by spatial configurations of attractants and repellents, and results of 
computation by structures of protoplasmic network formed by the plasmodium on 
the data sites [571 UHl IE] ■ The problems solved by plasmodium of P. polycephalum 
include shortest path [571 [33] , implementation of storage modification machines [S] , 
Voronoi diagram [37|, Delaunay triangulation logical computing [45l [6], and 
process algebra [3S]; see an overview in [B]. Previously [^ we have evaluated the 
road-modelling potential of P. polycephalum, however no conclusive results were 
presented back in 2007. A step towards biological approximation, or evaluation, of 
man-made road networks was done in our previous papers on the approximation 
of motorways/highways in the United Kingdom f9^, Mexico [T0| and Australia [12] 



Date: December 21, 2011. 



2 



ADAMATZKY, DE BAETS, AND VAN DESSEL 



by Plasmodium of P. polycephalum. For these countries we found that, in princi- 
ple, the network of protoplasmic tubes developed by plasmodium matches, at least 
partly, the network of man-made transport networks. However a country's shape 
and spatial configuration of urban areas, which are experimentally represented by 
sources of nutrients, might play a key role in determining the exact structure of 
Plasmodium networks. Also we suspect that the degree of matching between the 
Physarum networks and the motorway networks is determined by original govern- 
ment designs of motorways in any particular country. This is why it is so important 
to collect data on the development of plasmodium networks in all major countries, 
and then undertake a comparative analysis. 

Belgium is a good testbed for the evaluation of slime-mould approximation of 
motorways because 

• Belgium is an artificial country created relatively recently, in 1830. 

• It is amongst the most populated area in Europe. 

• There is a density misbalance between two major communities: Flanders 
is more densely populated than Wallonia. 

• the Belgian economy is centred around Brussels, by far biggest city, with 
hundreds of thousands of workers commuting to Brussels every day. 

In the early days the Belgian highways were constructed to provide a solution 
against the overcharged national and local roads, caused by the expanding number 
of cars. In the North of the country, construction was generally based on grow- 
ing demands from the economic and touristic sectors. The first highway was the 
one between Brussels and Ostend. Another 'early' highway was the one between 
Antwerp and Liege (E313) to open up the port of Antwerp's access to the 'hinter- 
land'. At the end of 1972 the most important cities were connected by highways. 
However, from the point of view of transport economics, only two of them were 
answering to an economic demand, a need for construction based on increasingly 
busy roads: the one between Brussels and Antwerp (E19), and the one between 
Brussels and Liege (E40). The others were intended as an investment trigger. The 
Autoroute de Wallonie (E42) was aimed at the economic reconversion of the old 
industrial axis in Wallonia (steel and coal industry). The E17 and E34 motorways 
provided an additional connection between the port of Antwerp and the French 
and German inner lands. The purpose of the E314 was to open up the province of 
Limburg (coal mining industry), and to provide a shortcut between Antwerp and 
the German Rhineland (Ruhrgebiet). Highway construction has been the result of 
political negotiations and the desire or need of the Northern and Southern partners 
to balance large investments in both parts of the country [51] . 

The paper is structured as follows. In Sect. [2] we give an overview of the ex- 
perimental setup employed. Analysis of protoplasmic networks produced by slime 
mould P. polycephalum in laboratory experiments is provided in Sect. [3] We com- 
pare slime mould generated and man-built motorways in Sect. [4] SectionjS] considers 
protoplasmic networks and Belgian motorways in the context of planar proximity 
graphs. Relations between the experimental results and the administrative subdivi- 
sion of Belgium are discussed in Sect. |6] In Sect. [7] we study outcomes of large-scale 
contamination and resulting reconfiguration of slime mould transport networks. 
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2. Experimental 

Plasmodium of P. polycephalum is cultivated in plastic container, on paper 
kitchen towels moistened with still water, and fed with oat flakes. For experiments 
we use 120 x 120 mm polystyrene square Petri dishes and 2% agar gel (Select agar, 
by Sigma Aldrich) as a substrate. Agar plates, about 2-3 mm in depth, are cut in 
the shape of Belgium. 

We consider the twenty-one mosl[^ populous urban areas in Belgium U (Fig. [T^), 
shown below in descending order of population size: 



(1) 


Brussels area, including Dilbeek 


(11) 


Mechelen 




and Vilvoorde 


(12) 


Kortrijk area, including 


(2) 


Antwerp area, including Beveren 




Mouscron and Waregem 




and Brasschaat 


(13) 


Hasselt 


(3) 


Gent 


(14) 


Oostende 


(4) 


Charleroi area, including La 


(15) 


Sint-Niklaas 




Louviere and Chatelet 


(16) 


Tournai 


(5) 


Liege area, including Seraing, 


(17) 


Genk area, including Maas- 




Verviers and Herstal 




mechelen 


(6) 


Brugge 


(18) 


Roeselare 


(7) 


Namur 


(19) 


Turnhout 


(8) 


Leuven 


(20) 


Arlon 


(9) 


Mons 


(21) 


Sankt-Vith 


(10) 


Aalst 







To represent areas of U we place oat flakes in the positions of agar plate corre- 
sponding to the areas. At the beginning of each experiment an oat flake colonised 
by Plasmodium is placed in the Brussels area (Fig. [2]). Our choice of inoculation site 
does not reflect the historical development of transport routes in Belgium (where 
inoculation should start in Aalst, Brugge, Kortrijk, or Gent) however it conveys 
the overwhelming economic power of the capital. We undertook 28 experiments. 
The Petri dishes with plasmodium are kept in darkness, at temperature 22-25°C, 
except for observation and image recording. Periodically, the dishes are scanned 
with an Epson Perfection 4490 scanner. 

3. SlIME MOULD TRANSPORT NETWORKS: BIO-ESSENTIAL MOTORWAYS 

Plasmodium is inoculated in the Brussels area. In the first 24 h it propagates 
towards and occupies Leuven and Mechelen, and then propagates from Mechelen 
to Antwerp, from Antwerp to Sint-Niklaas (Fig. [3^b). In the next 24 h slime 
mould propagates from Sint-Niklaas to Aalst and Gent, from Gent to Brugge and 
from Aalst to Kortrijk. Links from Brugge to Oostende and Roeselare are built 
during the same time interval (Fig. [3];d). Westward development of plasmodium is 
somehow stopped. Despite an attempted propagation from Turnhout towards the 
Hasselt and Genk areas the slime mould never actually reaches these areas in the 
first 48 h (Fig.|3j:d). 

By the 72nd hour after being inoculated in Brussels almost all urban areas but 
Hasselt and Genk are colonised by slime mould. Namely in the time interval 48-72 h 



Arlon and Sankt-Vith are not amongst the most populated areas but we added them for 
completeness, to allows the slime mould propagating towards Luxembourg and Germany 
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Figure 1. Experimental setup: (a) outline map of Belgium [21] 
with major urban areas U shown by encircled numbers; (b) urban 
areas, represented by oat flakes, are colonised by slime mould of 
P. polycephalum. 
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Figure 2. Oat flakes, representing urban areas U, colonised by 
slime mould. The growing substrate is on top of the Belgian pop- 
ulation density map |39j . 



ADAMATZKY, DE BAETS, AND VAN DESSEL 




(c) 72 h (f) 72 h 



Figure 3. Example of anti-clockwise colonisation, where south 
of Belgium is colonised via Tournai— >-Mons— >-Charleroi— >-Namur 
route. (ace) colour scans of the experimental Petri dishes, 
(bdf) colour enhanced and grey-scaled versions. Time lapsed after 
slime mould inoculation in Brussels is shown in captions to sub- 
figures. 



SLIME MOULD IMITATION OF BELGIAN MOTORWAYS 



7 




(a) 96 h (b) 96 h 

Figure 4. Final stages of anti-clockwise colonisation Fig.|3] 



Plasmodium grows from Aalst to Tournai, from Tournai to Mons, and from Mons 
to Charleroi. Slime mould branches at Charleroi and grows in parallel to Namur 
and Arlon. It propagates from Namur to Sankt-Vith, from Sankt-Vith to Liege. 
(Fig. |3]3f). By the 96th h after being inoculated the slime mould propagates from 
Liege to Genk and Hasselt areas, and from Hasselt to Leuven. The Plasmodium's 
explorative activities in the Brussels area are 'resumed' when the slime propagates 
from Leuven and re-occupies the Brussels area (Fig. |4] and Fig. [§1.). 

Colonisation of Belgium by slime mould shown in Fig. [5] develops initially ac- 
cording to the scenario described above with the following deviations. Plasmodium 
propagates from Gent to Brugge, Roeselare and Kortrijk at the same time (Figs. [5^- 
d andjoja). Slime mould grows from Leuven to Namur and then to Charleroi. From 
Charleroi plasmodium propagates to Mons and from Mons to Tournai (Figs.[5^b). 
Arlon and Sinkt-Vith are reached by slime mould via Leuven, Hasselt and Liege 
(Figs. [5]:d and[5^f). Antwerp is never colonised by slime mould in this particular 
experiment (Fig. [6]). 

Two examples discussed above show that the dynamics of colonisation varies 
from experiment to experiment. Therefore from experimental data we extract a gen- 
eralised Physarum graph. To generalise our experimental results we constructed a 
Physarum graph with weighted edges. A Physarum graph is a tuple P = (U, E, w), 
where U is a set of urban areas, E is a set of edges, and w : E — >• [0, 1] associates 
to each edge of E a probability (or weights) . For every two regions a and b from U 
there is an edge connecting a and 6 if a Plasmodium's protoplasmic link is recorded 
at least in one of k experiments, and the edge (a, b) has a probability calculated as 
the ratio of experiments where protoplasmic link (a, b) occurred in the total number 
of experiments (fc = 23). For example, if we observed a protoplasmic tube connect- 
ing areas a and 6 in 7 experiments, the weight of edge {a,b) will be w{a^b) = 
We do not take into account the exact configuration of the protoplasmic tubes but 
merely their existence. Further we will be dealing with threshold Physarum graphs 
P(6') = (U, r(E), w, 0). A threshold Physarum graph is obtained from Physarum 
graph by the transformation: r(E) = {e G E : w{e)geq9}. That is, all edges with 
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(c) 96 h (f) 96 h 



Figure 5 . Example of clockwise colonisation of south of Belgium 
with dominating route Leuven— >-Hasselt— >-Liege. (ace) colour scans 
of the experimental Petri dishes, (bdf) colour enhanced and grey- 
scaled versions. 
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Figure 6. Scheme of colonisation dynamics illustrated in 
(a) Fig. |3] and (b) Fig. [5| 



weight less than are removed. Examples of threshold Physarum graphs for various 
values of are shown in Fig. [THj 

The 'raw' Physarum graph P(^) is a non-planai|^ acyclic graph (Fig. [t^l). In 
a raw graph each edge appears at least in one laboratory experiment. We call an 
edge of Physarum credible if this edge is represented by a protoplasmic tube in over 
20% of laboratory experiments. Such graph P(^) is shown in Fig. [Tja. It is still 
non-planar, however, the only intersecting edges are the links Liege area — Arlon 
and Namur — Sankt-Vith. 

Finding 1. The Physarum graph with all credible edges is a non-planar cyclic 
graph. 

The graph remains connected while 9 grows up to || (Fig. . The Turnhout 
urban area becomes an isolated vertex when = ^ (Fig- ■ For this value of 
9 the Physarum graph becomes planar because the link (Liege area - Arlon) is 
represented in more experiments than the link (Namur - Sankt-Vith). The graph 
P(||) has the largest (among all Physarum graphs studied here) empty, i.e. not 
having any edges inside, circle. Clockwise, starting in the Brussels area it spans 
Leuven, Hasselt, Liege, Namur, Charleroi, Mons, Tournai, Kortrijk, Gent, Aalst 
and finishes in Brussels. 

The Physarum graph P(6') splits into three components when 9 = The 
smallest component is the isolated vertex Turnhout. The medium size component 
is a cycle Liege — Hasselt — Genk — Liege attached to a segment Arlon — Sank- 
Vith (Fig. [T^). The largest component is a proximity graph spanning all remaining 
urban areas. 

Finding 2. The following slime mould transport chains appear in almost all ex- 
periments (Fig.^): 

• chain Roeselare — Kortrijk — Mons — Charleroi — Namur 



A planar graph consists of nodes which are points of the Euchdcan plane and edges which are 
straight segments connecting the points. 




(c) e = i (d) e = i| 




(e)e = i {f)e = i 

Figure 7. Physarum graphs P(^) for critical values of 6. 

• chain Oostende — Brugge — Gent — Aalst — Brussels — Leuven — 
Mechelen -— Antwerp Sint-Niklaas 

• chain Liege — Sankt-Vith — Arlon 

• one-link chain Hasselt — Genk. 
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4. Slime mould versus Belgian motorways 

As we can see in the examples of the experimental configurations in Fig. [8j 
Plasmodium networks are polymorphic and no two networks are exactly the same. 
Some motorways are matched by protoplasmic tubes well, others just approximated, 
and some do not have a slime mould representation at all. 

For example. Fig. [8^ demonstrates that slime mould developed protoplasmic 
tubes corresponding to motorway Af 0/E40 between Gent and Brugge, Ef 7 between 
Gent and Kortrijk, E17 and E403/Af7 between Gent and Tournai, and E42 con- 
necting Mons to Charleroi to Namur to Liege. Motorways E42 (Tournai - Mons), 
E46 (Antwerp - Turnhout) and E40 (Leuven - Liege) are represented by slime 
mould in the experiment illustrated in Fig. [8]d. At the same time, the transport 
link E4f f from Brussels to Namur to Arlon is not represented in the configurations 
inFig.|8} 

Due to the variability of the slime mould networks it would be unreasonable to 
attempt comparing the exact topology of slime mould and man-made transport 
links. Instead we will compare them at the level of graph-theoretical representa- 
tions. 

The graph H of the Belgian motorway network (Fig. [9]) is constructed as follows. 
Let U be a set of urban regions/cities; for any two regions a and b from U, the 
nodes a and b are connected by an edge (a, b) if there is a motorway starting in 
the vicinity of a, passing in the vicinity of 6, and not passing in the vicinity of any 
other urban area c € U. In the case of branching - that is, a motorway starts in a, 
goes in the direction of b and c, and at some point branches towards b and c - we 
then add two separate edges (a, b) and (a, c) to the graph H. The motorway graph 
is planar (Fig. |9]) . 

Finding 3. Physarum polycephalum almost completely approximates the Belgian 
motorway network. 

Namely, H ^ P(^) but 25 of 28 edges of H are edges of P(^). The following 
motorway links are never represented by protoplasmic tubes, not in a single exper- 



iment: Brussels to Tournai, Brussels to Antwerp, Antwerp to Hasselt (Fig. 10 a, and 
Fig.[9|. 

Finding 4. On redundancy. Motorway links connecting Brussels with Antwerp, 
Tournai, Mons, Charleroi, Namur, and links connecting Leuven with Liege and 
Antwerp with Genk and Turnhout are proved to be redundant components of the 
Belgian transport system in the slime mould experiments. 

The skeletal motorway network represented by the plasmodium network in lab- 
oratory experiments is shown in Fig. |10^ 3. This network is represented by pro- 
toplasmic tubes in over 35% of experiments; and in almost 40% of experiments 
without the Antwerp — Turnhout link (Fig. [Toj^d) . By increasing the level of slime 
mould's "confidence" to almost 60% we loose motorway links connecting Leuven 



with Hasselt, Namur with Liege and Arlon, and Liege with Arlon (Fig. 10 ;d) 



Finding 5. Motorway segments All (Antwerp — Sint-Niklaas), A19 (Antwerp 

— Mechelen), E42 (Liege — Sankt-Vith), AIO/E4O (Oostende — Gent — Aalst 

— Brussels — Leuven), A17/E403 (Roeselare — Kortrijk — Tournai), E42/E19 
(Tournai — Mons), E42 (Mons — Charleroi — Namur) are essential transport 
links from the slime mould's point of view. 




Figure 8. Examples of P. polycephalum protoplasmic networks 
on the Belgian motorway network. 
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Figure 9. Graph H of Belgian motorway network. 



As illustrated in Fig. [TOf these transport links are represented by the protoplas- 
mic network in almost 80% of laboratory experiments. 

5. Proximity graphs 

A planar proximity graph is a planar graph where two points are connected by 
an edge if they are close in some sense. A pair of points is assigned a certain neigh- 
bourhood, and points of the pair are connected by an edge if their neighbourhood 
is empty. Here we consider the most common proximity graph as follows. 

• GG: Points a and h are connected by an edge in the Gabriel Graph GG 
if the disc with diameter dist{a, b) centred in middle of the segment ab is 
empty [BIIS] (Fig. [5^). 

• RNG: Points a and b are connected by an edge in the Relative Neighbour- 
hood Graph RNG if no other point c is closer to a and b than dist{a, b) |31] 
(Fig#). 

• MST: The Euclidean minimum spanning tree (MST) |30j is a connected 
acyclic graph which has the minimum possible sum of edges' lengths (Fig.jsjD) 

In general, the graphs relate as MST C RNG C GG [201 [551 gl] ; this is called the 
Toussaint hierarchy. 

Finding 6. The Belgian motorway graph is approximated by the Gabriel graph with 
80% accuracy and by the relative neighbourhood graph with 70% accuracy. 

The following motorways links from H are not presented in the Gabriel graph 
GG: (Brussels area — Antwerp area), ( Brussels area — Tournai), (Brussels area 
— Namur), (Antwerp area — Hasselt), (Leuven — Liege area), (Liege area — 
Arlon) (Fig. [T2|i). With regards to the relative neighbourhood graph, a few more 
transport links from H are not a part of RNG: (Brussels area — Mons), (Brussels 



14 ADAMATZKY, DE BAETS, AND VAN DESSEL 




(e) HnP(i) 

Figure 10. Intersection H n P(0) of motorway graph H with 
Physarum graphs P(6') for critical values of 9. 



area — Leuven) and (Namur — Arlon) (Fig. [T2]d). Thus GG represents 23 of 29 
edges of the motorway graph H and RNG 20 of 29 edges of H. 

Finding 7. The minimum spanning tree rooted in Brussels is not a subgraph of 
the Belgian motorways graph, MST ^ H. 
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(a) RNG 



(b) GG 




(c) MST 

Figure 11. Proximity graphs constructed on sites of U. (a) rela- 
tive neighbourhood graph, (b) Gabriel graph, (c) minimum span- 
ning tree rooted in Brussels. 



This is because the links (Mechelen — Sint-Niklaas) and (Mechelen — Leuven) 
exist in MST but not are not present in H (Fig. [l2);). The minimum spanning 
tree is considered to be an optimal, in a sense of minimality of edge lengths, acyclic 
planar graph. The fact that two of the spanning tree edges are not represented by 
man-made motorway links allows us to suggest that the Belgian motorway network 
is not optimal, at least not optimal in spanning of major urban areas. 

Considering Physarum graph P(|g) is credible because its edges are represented 
by protoplasmic tubes in at least 40% of laboratory experiments, we can compare 
it with three basic proximity graphs (Fig. 13). 



Finding 8. If slime mould represented A7/A54 Brussels — Charleroi, non-existing 
in reality motorways between Leuven — Namur and Turnhout — Hasselt its credible 
transport network would be a super-graph of the relative neighbourhood graph. 

This is a direct outcome of comparing Fig. [T3^ and Fig. [TT^ . Relative neigh- 
bourhood graphs are considered to be optimal cyclic graphs in terms of total edge 
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(a) RNG n H 



(b) GG n H 




(c) MST(l) nH 

Figure 12. Intersections of the Belgian motorway graph with 
(a) relative neighbourhood graph, (b) Gabriel graph, (c) minimum 
spanning tree rooted in Brussels. 



length and travel distance, and are known to be a good approximation of road 
networks [321 HOj • The Leuven — Namur and Brussels — Charleroi links are rep- 
resented by slime mould in less than 20% of laboratory experiments (Fig.[7|D). The 
intersection of the Physarum graph P(||) (Fig. [7|i) with the minimum spanning 
tree (Fig. 11;) comprises two disconnected components: one lies in Flanders and 



another in Wallonia (Fig. 13' 



6. Dissolution: Transport networks and administrative subdivision 

The province Brabant Wallone has no transport routes originating in and passing 
through it for 6* = jf- When increases to the Antwerp province becomes 
isolated from the other provinces of Belgium. For = || we have isolated Antwerp 
and Limburg provinces, and clusters of interconnected regions: 

(1) West-Vlaanderen, Hainaut and Namur provinces, 

(2) West-Vlaanderen, Oost-Vlaanderen, Vlaams Brabant, Oost-Vlaanderen provinces. 
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(a) RNGnP(§) 



(b) GGnP(§) 




(c) MSTnP(ll) 



Figure 13. Intersections of the Physarum graph P(5|) with 
(a) relative neighbourhood graph, (b) Gabriel graph, (c) minimum 
spanning tree rooted in Brussels. 



(3) Liege and Luxembourg provinces. 



In terms of spanning tree (Fig. 11 ) the transport link between the Brussels area 
and the Charleroi area is the only means of keeping the Dutch and French speaking 
communities together. 



Finding 9. // the two parts of Belgium were separated with Brussels in Flanders, 
the Walloon region of the Belgian transport network would be represented by a single 
chain from Tournai in the north-west to the Liege area in the north-east and down 
to southmost Arlon. 



This transport link is not part of the spanning tree rooted in Brussels, and when 
omitted the two Belgian communities become isolated from each other (Fig. 14 1. 
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Flanders 




Figure 14. Intersection of Physarum graph P(||) with the min- 



imum spanning tree rooted in Brussels. 



7. Response of slime mould networks to imitated disasters 

To study the reaction of Physarum-grown transport networks on major disasters, 
we placed crystals of sodium chloride in the approximate positions of Doel Nuclear 
Power Plant, near Antwerpen (seven experiments) and Tihange Nuclear Power 
Station, near Liege (nine experiments). Sodium chloride is a strong repellent for P. 
polycephalum. The salt diffuses in the substrate outwards its original application 
site (epicentre of disaster). It can therefore imitate radioactive and/or chemical 
pollution and subsequent disturbance spreading along Belgian transport networks. 
Images of protoplasmic networks reconfigured 24 h after start of contamination are 
shown in Figs. [17] and [18] 

A typical response to propagating contamination is dissected in Fig. |l5] In this 
example we place a salt grain on the oat flake representing the Liege area (Fig. 15 
'a'). In 24 h after the imitated accident, the contamination spreads as far Namur on 
the west and Sankt-Vith on the east and Hasselt and Genk in the south. Transport 
links in proximity of the contamination epicentre become destroyed and abandoned. 



An example of abandoned transport link is a protoplasmic tube (marked 'b' Fig. 15 ) 



representing E42/A27 (Liege to Sankt-Vikt motorway). No alternative routes for 
such destroyed transport links are offered by the slime mould. 

Transport links being at a significant distance from the epicentre but yet inside 
the contamination zone are shifted further away from the contamination. For exam- 
ple, slime mould abandons protoplasmic tube representing motorway E411, marked 
'c' in Fig. 15 but grows another tube, marked 'd' in Fig. 15 slightly westwardly. 



Urban areas not directly affected by contamination show signs of increased explo- 
rative and scouting activity, for example, Antwerp and Turnhout areas, marked 
'e' in Fig. [15] Also transport links not affected by contamination became visibly 
enhanced (Fig. 15 'g'). In the situation of continuing contamination plasmodium 



'considers' the last opportunity to survive by forming sclerotium (hardened body of 
the 'hibernating' slime mould) , initial stage of the sclerotium formation is labelled 
'f ' in Fig. [151 



Based on outcomes of our scoping experiments (Figs. 17 and 18 ) we can propose 
the following scenarios of response to contaminations: 
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Figure 15. Illustration of the slime mould's response to contam- 
ination, (a) epicentre of contamination, (b) abandoned and not 
repositioned transport link, (c) abandoned and repositioned trans- 
port link, (d) new location of the abandoned transport link, (e) ex- 
plorative activity, (f) emergency preparations, initial stage of scle- 
rotinisation. (g) enhanced transport links. Image is taken 24 h 
after initiation of contamination in Liege area. 



Epicentre of contamination is at the Tihange Nuclear Power Station near 
Liege. Segments of transport network connecting Namur, Liege, Hasselt, 
Genk, Sankt-Vith and Arlon are destroyed or abandoned. Hyperactivity 
is observed in domains surrounding Sinkt-Niklaas, Antwerp and Turnhout. 
Preparations for emergency hibernation take place at Aalst, Brussels, Leu- 
ven and Mechelen. Transport links between Gent, Roeselare, Kortijk, Tour- 
nai, Mons, Charleroi, Namur are hypertrophied (Fig. 18). 
Epicentre of contamination is at the Doel Nuclear Power Plant near Antwer- 
pen. Domains surrounding Antwerp, Sinklaas, Mechelen, Turnhout become 
depopulated and transport links are abandoned. Explorative activity is ob- 
served in areas of Oostende, Brugge, Roeselare, and Arlon and Sankt-Vith. 
Attempted mass-migration is recorded into Northern France, South of the 
Netherlands and Luxemburg and West Germany. Transport links in the 
chain Brugge — Roeselare — Kortijk — Tournai — Mons — Charleroi 
— Namur — Liege — Sankt-Vith/ Arlon are significantly hypertrophied 



(Fig. 17 1. Preparations for emergency hibernation take place at Brussels 



and Leuven, and at Hasselt and Genk areas. 



20 



ADAMATZKY, DE BAETS, AND VAN DESSEL 



8. Conclusions 

To evaluate how good Belgian motorways are from an amorphous living creature 
point of view, we conducted the following laboratory experiments with slime mould 
P. polycephalum. We represented major urban areas with oat flakes and inoculated 
slime mould in the oat flake corresponding to Brussels. We waited till the slime 
mould colonised all oat flakes and then analysed the slime mould's protoplasmic 
network structure and compared it with the man-made motorway network and 
basic planar proximity graphs. We found that P. polycephalum almost but not 
completely approximates the Belgian motorway network. Transport links Roeselare 

— Kortrijk — Mons — Charleroi — Namur, Oostende — Brugge — Gent — Aalst 

— Brussels — Leuven — Mechelen — Antwerp — Sint-Niklaas, Liege — Sankt-Vith 

— Arlon, and Hasselt — Genk appear as protoplasmic tubes in almost all laboratory 
experiments. Motorway links connecting Brussels with Antwerp, Tournai, Mons, 
Charleroi, Namur, and links connecting Leuven with Liege and Antwerp with Genk 
and Turnhout are "considered" by slime mould to be redundant and thus almost 
never appear in our experiments with the slime mould. If slime mould represented 
A7/A54 Brussels — Charleroi, non-existing in reality, motorways between Leuven 

— Namur and Turnhout — Hasselt its credible transport network would be a super- 
graph of the relative neighbourhood graph. 

Motorway segments A17 (Antwerp — Sint-Niklaas), A19 (Antwerp — Meche- 
len), E42 (Liege — Sankt-Vith), AI0/E40 (Oostende — Gent — Aalst — Brussels 

— Leuven), A17/E403 (Roeselare — Kortrijk — Tournai), E42/E19 (Tournai — 
Mons), E42 (Mons — Charleroi — Namur) are essential transport links from the 
slime mould's point of view. If the two parts of Belgium were separated with Brus- 
sels in Flanders, the Walloon region of the Belgian transport network would be 
represented by a single chain from Tournai in the north-west to Liege area in the 
north-east and down to southmost Arlon. 

While imitating major disasters leading to contamination propagating from Liege 
and Antwerp areas we found several scenarios of major restructuring of transport 
networks and possible routes of mass-migration. 

We believe the results of our scoping experiments and analysis can be used not 
only in bio-inspired unconventional computing but also in more classical fields of 
urban and transport planning. Possible applications include but are not limited to 
restructuring of rural landscape 17J, novel approaches towards mapping vehicular 
accessibility '46] , simulation of traffic dynamics in Belgian motorways |14| , develop- 
ment of Trans-European transport networks |15j , bottom up landscape planning |35j 
and modelling relations between urban sprawl and growing transport networks |32j . 
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9. Appendix A: Physarum graphs for all values of 9 
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Figure 17. Response to spreading contamination recorded in lab- 
oratory experiments. Epicentre of contamination is in Antwerp 
area. 
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Figure 18. Response to spreading contamination recorded in lab- 
oratory experiments. Epicentre of contamination is in Liege area. 



